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Aims. A study of the six largest coronal X-ray flares in the Chandra archive is presented. The flares were observed on II Peg, OU And, Algol, 
HR 1099, TZ CrB and CC Eri, all with the High Energy Transmission Grating spectrometer (HETG) and the ACIS detectors. 
Methods. We reconstruct an Emission Measure Distribution EMD(T), using a spectral line analysis method, for flare and quiescence states 
separately and compare the two. Subsequently, elemental abundances are obtained from the EMD. 

Results. We find similar behaviour of the EMD in all flares, namely a large high-r component appears while the low-T (kT < 2 keV) plasma 
is mostly unaffected, except for a small rise in the low-J Emission Measure. In five of the six flares we detect a First Ionization Potential (FIP) 
effect in the flare abundances relative to quiescence. This may contradict previous suggestions that flares are the cause of an inverse FIP effect 
in highly active coronae. 
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1. Introduction 

The ongoing multi-wavelength research on stellar coronae has 
shown that coronal activity in cool stars is closely related to 
the magnetic activity. While in the solar case the coronal struc- 
tures are spatially resolved, one has to rely on accurate spec- 
troscopy in analyzing the unresolved stellar corona. Although 
showing general similarity to the sun, some stellar systems ex- 
hibit coronal activity 3-4 orders of magnitude more energetic 
than in the sun, especially during large flares. Whether this in- 
dicates a scaled up version of the solar coronal activity or a 
differently structured corona is not yet clear. 

Another effect which is not well understood is the abun- 
dance variations between the photospheres and coronae of 
the stars. In the solar case this has been known as the First 



possibly by flares that evaporate high-FIP material from the 



lower chromosphere to the corona (Brink man et all l2001) and 



by electric currents t hat suppress the diff usion of low-FIP 



species to the corona (ITelleschi et al.l 2005). Observationally, 



however, the associatio n to date o f flares with FIP or IFIP 
effects is ambiguous. iGiidel etalJ dl999h and Audard et al 



(2001) found an increase in low FIP abundances during 
flares on the RSCVn UX Ari and HR 1099, respectively. In 
some cases the variations in abundances were not FlP-related 



Ionization Potential (FIP) effect (Feldman, 19921) . where com- 
pared to their photospheric abundances, elements with low FIP 
(under~10 eV) are over-abundant relative to elements with 
higher FIP. Many systems show this effect, but some active 
systems show an opposite effect where it is the high FIP el- 
ements who are over abundant. This was labeled the Inverse 



FIP (IFIP) effect (Brin kman et al. . 2001). Audard et al. (2003) 



found that highly active RS CVn binaries show an IFIP effect 
while less active systems sh ow either no effect or a solar FIP 
effect. ITelleschi et al.l y005) found a related result in a sam- 
ple of solar like stars, where abundances change from IFIP to 
FIP with the age (and decreasing activity) of the star. This has 
led to the suggestion that activity affects coronal abundances, 



(Ost en et aU l2003; Gtid el et aUl2004l) . In other cases, n o abun- 
dance variations were found at all, e.g., Mag gio et all (bOOOh 
and iFranciosini et al.l ((2001), though the last two works are 
based on low spectral resolution data. 

The X-ray band offers several advantages in the study of 
stellar flares. The high temperatures in the flares, typically with 
kT > 1 keV, mean that most of the emission is in the X-ray 
band and the major emission lines of the highly ionized ele- 
ments fall in the 1 .8 to 40 A range. This provides measurements 
of line fluxes from up to 10 Fe K- and L- shell ionization de- 
grees (Fe XVII to Fe XXVI), as well as from the two K-shell 
ions of the other common elements. While previous instru- 
ments limited plasma emission models to two- or three- thermal 
components, line resolved spectra from Chandra enable the re- 
construction of a more accurate distribution of plasma temper- 
atures, as well as better measurements of electron densities and 
abundances. 

The reconstruction of the emission measure distribution 
{EMD) that describes the plasma thermal structure, and the 
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Table 1. Observations used in this work. 



Obs. ID 


HD 


Other Name 


Exposure (ks) 


Start Date 


Type 


Radius 


Distance (pc) H 


Pari, (d) 


1451 


224085 


II Peg 


43.3 


1999-10-17 23:28:28 


K2IV+? 


2.2 5 


42.34 


6.7 


1892 


223460 


OU And 


96.9 


2001-08-11 00:18:00 


Single GUIIe 


9 F 


135 




604 


19356 


Algol 


52.4 


2000-04-01 02:20:34 


B8V+G8III S 


2.8/3.54* 


28.46 


2.86 B 


62538 


22468 


HR1099 


95.9 


1999-09-14 22:53:10 


G5IV+K1IV S 


1.3/3.9 5 


28.96 


2.84 


15 


146361 


TZCrB 


84.8 


2000-06-18 13:41:55 


F6V+G0V S 


1.22/1.21* 


21.69 


1.14 


6132 


16157 


CCEri 


30.95 


2004-10-01 01:46:49 


K7Ve/M4 s 


ojp s 


11.51 


1.56 


4513 


16157 


CCEri 


89.45 


2004-10-01 20:54:39 











iPerrymanet al,~1997) 



H 

S - 

B 

F 



The HIPPARCOS catal og 
Strassmeier et al.| Jl993l) 
Budding etalJd2004h 



Fekel et al.Ul986l) 



measurements of abundances are difficult, not only because 
the two are entangled, but also due t o the poor mathematica l 
definition of the problem as shown bv lCraig & Brown J1976I) . 
Small variations of observed line fluxes result in large varia- 
tions in the derived EMD. Almost any variation of the EMD 
can be accommodated to produce the same spectra, up to mea- 
surements uncertainties, if done on small enough temperature 
scales. Consequently, many authors avoid setting confidence 
intervals on their fitted EMD, or use EMD smoothing, which 
makes the comparison of EMDs produced by different meth- 
ods a difficult task. The correlation between EMD parameters 
and deduced abundances is often neglected or ignored. 

In this work, we scanned the public Chandra archive for 
bright flares that enable good line flux measurements. Six such 
flares were found on the systems: II Peg, OU And, Algol, HR 
1099, TZ CrB & CC Eri. While these are not the only flares 
in the archive, they represent the brightest and best resolved 
flares available. Some of these observations have been analyzed 
before, but not always with regard to the flare and not in a way 
that allows a comparison with other analysis methods. We refer 
to these previous works in section 15721 

We apply the same analysis methods, employing the deriva- 
tion of a continuous EMD and relative abundances, for the six 
targets in flaring a nd quiescence states . The method is based 



on the one used in lNordon et a l. (2006). Emphasis is given to 



the ability to compare the results in a statistically significant 
way in order to seek real variations in the thermal structure and 
abundances between flare and quiescent states. In particular, 
we want to establish whether the large flares had a statistically 
significant effect on the coronal abundances and whether it is 
FlP-related. 

2. Observations 

2.1. The Sample 

The Chandra public archive was searched for observations of 
cool stars (spectral types A to M) in any system configuration, 
featuring strong, long-duration flares. The main selection cri- 
terion was the requirement for enough photons in the flare to 
allow a detailed spectral analysis. So, while there are other 
observed flares not included in this work, this sample repre- 
sents the best flares, in terms of photon counts, observed with 
Chandra. All grating observations were examined, however, 



the selected sample happens to contain only the ACIS+HETG 
instrument configuration, which is the most common configu- 
ration used for coronal targets. 

The details of the selected observations and targets are sum- 
marized in table Q] The data for CC Eri are composed of two 
observations separated by a short time gap. The flare occurred 
at the end of the first observation and the last part of the flare 
decay was cut-off. Due to low flux in the quiescence state, we 
integrated both the period before the flare from the first obser- 
vation, and the entire second observation for the spectral ex- 
traction of the quiescence state. 

2.2. Light Curves and Spectra 

Light curves were produced using combined counts of the High 
(HEG) and Medium (MEG) energy grating arms of all orders as 
well as the zero-order region. Background, though negligible, 
was estimated and subtracted using off-source CCD regions. 
The light curves are presented in Figure Q] using 400 s bins. 
Segments used for flare and quiescence spectra extraction are 
marked with F and Q respectively. 

Flare and quiescent spectra are presented in figures|2]-[3]and 
correspond to the time segments of the observations marked in 
figure Q] In the Algol observation, the time segment before the 
flare was excluded as the system was still in eclipse. The plots 
are in 0.01 A bins and use combined fluxed spectra of HEG and 
MEG. 

3. Modelling Method 

3.1. Emission measure distribution 

While seeking the plasma parameters that reproduce the prop- 
erly measured fluxes of several selected lines, we are particu- 
larly interested in examining the thermal structure and elemen- 
tal abundances independently. Therefore, we develop and use 
a method that disentangles the mutual dependence in a simple 
way. The continuum emission is problematic for two reasons; 
One being the typical low level of continuum emission in quies- 
cent states. The other reason is the lack of unique features in the 
continuum, which makes the distinction between thermal and 
non-thermal components, or between high-T components and 
low metalicity, difficult and ambiguous. Better constraints are 
provided by the high energy end of the bremsstrahlung spec- 




Time Observing [ks] Time Observing [ks] 



Fig. 1. Light curves of the observations using 400 s bins. All orders of diffraction are included and background is subtracted. 
The vertical dashed lines indicate the segments considered as flaring or quiescence and marked with letters F and Q respectively. 
The CC Eri observation consists of two nearly continuous observations. 



trum, but in cases of a wide temperature distribution even the 
distinctive bremsstrahlung peak may be blurred. Moreover, in 
flares, the continuum turnover is often beyond the instrument 



band. Thus, we avoid using the continuum and rely on the more 
accurate line fluxes. 

The observed line flux F q . of ion q due to the atomic tran- 
sition j — > i can be expressed by means of the element abun- 
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Fig. 2. Time averaged flare and quiescence spectra of the 6 targets in the 1.7-12 A range. HEG and MEG gratings 1st orders are 
combined at MEG resolution and 0.01 A bins. Orange is flare and blue is quiescence. 
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Fig. 3. Time averaged flare and quiescence spectra of the 6 targets in the 12-20 A range. HEG and MEG gratings 1st orders are 
combined at MEG resolution and 0.01 A bins. Orange is flare and blue is quiescence. 
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dance with respect to hydrogen A z , the distance to the object 
d, the line power P 9 .., the ion fractional abundance f q and the 
EMD as: 



n 



And 2 



P q ji(T)f q {T)EMD{T)dT 



(1) 



The line emissivity e*(T) = P 9 i (T)f q (T) describes ion and line 
specific parameters, while the general plasma parameters are 
described by the EMD as: 



EMD = n e n H dV/dT 



(2) 



Where n e and are the electron and hydrogen number den- 
sities averaged over the volume of plasma in the temperature 
interval [T, T + dT]. 

The strongest unblended line from each ion is selected 
for the analysis and density sensitive lines are avoided. 
Nonetheless, a complete set of ionic spectra are used to account 
for residual blending. The temperature dependence of the line 
emissivity eJ t (T) comes mainly from f q (T) and therefore a sec- 
ond line from the same ion contributes little additional infor- 
mation, except for more photon statistics. It could however, in 
some cases, provide statistical compensation for uncertainties 
in the atomic parameters of the strong line. 

We use the primary line from every Fe ion (Fe XVII to 
Fe XXV) in the observed spectra to get a set of integral equa- 
tions, whose solution yields the EMD scaled by the unknown 
Fe abundance: 



F q . = 



And 2 



f 

J T 



P q H {T)f q {T)EMD{T)dT 



(3) 



In order to include lines of other elements, but without the 
need to fit their abundances, we use ratios of the He-like to 
H-like line fluxes, from the same element. Thus, the element 
abundance A z cancels out. This adds another set of equations 
that constrain the shape of the EMD, but do not depend on the 
abundances: 



R- 



77 He- 1 ike 



pH-like 
t lk 



£~* Pf(T)f He (T)EMD{T)dT 
jT M P H (T)fH(T)EMD(T)dT 



(4) 



The X-ray spectra we use here include as many as ten Fe ions, 
but no more than two ions from other elements. For every ob- 
servation we are able to use similar equations, but different spe- 
cific ions, depending on which lines have a sufficiently good 
signal. 

The line fluxes and flux ratios are then fitted using the 
least squares best fit method to solve this set of equations 
for the EMD, where the EMD is expressed by a parameteric 
non-negative function of T (see details in section I3.31 l. The 
minimization algorithms used are Levenberg-Marquardt and 
Nedler-Mead (simplex) alternately and combined, taking the 
best result. The fit yields the estimated shape of the EMD, in- 
dependent of any assumptions for the abundances, and is scaled 
by the Fe abundance. In other words, we solve for the prod- 
uct of Af e ■ EMD. The integration in eqs. (|3]|I]l is started at 
kTo = 0.2 keV and cut-off at kT max = 10 keV, under and 



over which the EMD is completely degenerate, as no new lines 
emerge from these temperature regimes. This means that some 
of the flux in high-T lines may originate from plasma at even 
higher temperatures, the effect of which will be emission mea- 
sure (EM) added to the last high-T bin. The same applies to the 
lower temperature cut-off. Some emission in the low-T lines 
may originate from temperatures lower than the cut-off. For 
this reason, we avoid using lines with significant emissivity be- 
low the cut-off temperature, even if they are available (e.g. O 
VII). The coolest forming ion we use is O VIII whose emissiv- 
ity peaks at ~0.25 keV and f or which over 9 5% of the emissiv- 
ity lies above the threshold. |GueLaL (2006) showed that below 
0.2 keV HETG cannot constrain the EMD. The abundance- 
independent approach, namely using abundance independent 
equations for fitting the EMD and then using the EMD to re- 
solve the abun dances, has been u sed i n other works as well, for 
example: Ischmitt & Nesl d2004 and lTelleschi et all d2005l) . 

The atomic dat a for the line powers are calculated using 
the HULL AC code dBar-Shalom et alll200ll) . The ionic abun - 
dances (f„) for : Fe, A r, S, Si, Mg are taken from|Gu| J2003I) . 
Mazzotta et al.l (|1998) is used f° r the other elements. 



3.2. Line fluxes 

In order to measure the line fluxes and solve for possible line 
blending, we preform an ion-by-ion fitting to the spectra. The 
line powers for each ion are calculated at its maximum emis- 
sivity temperature and then passed through the instrument re- 
sponse. This accounts for all of the lines and fixes the ratios of 
lines from a given ion. 

The observed spectrum is then fitted by sets of com- 
plete individual-ion spectra simultaneously, together with a 
bremsstrahlung continuum, composed of several discrete- 
temperature components. This results in an excellent fit that 
accounts for all the observed lines and blends. Since we ulti- 
mately use only several selected lines for the analysis, we make 
sure that they fit well and the fitting of the other lines serves for 
estimating the blends in the selected lines (which is small to be- 
gin with). The fitted continuum parameters are independent of 
the reconstructed EMD and are not used in the following anal- 
ysis. The synthetic continuum serves only for measuring line 
fluxes above it. Due to the high resolving power of the HETG 
that produces narrow line profiles and the ionic spectra that in- 
clude all lines that may contribute to a pseudo continuum, the 
local continuum level is reliable. The narrow line profiles of 
HETG also mean that slight deviations in local continuum level 
have little effect on the measured li ne flux. This process is sim- 
ilar in principle to the one used in iBehar. et al.l (1200 lb and in 
Brinkman et alj d200ll) . 



The line fluxes measured from all the observations and used 
in the following analysis are listed in tables lA~TTlA.3l in the ap- 
pendix. 

3.3. EMD parametrization and fitting 

Our goal is to compare the EMD of the flare and quiescence 
states. It is important to note that the solution for the EMD 
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is n ot unique as is t he ca se with integral equations of this 
sort dCraig & Brown , 119761) . On scales much smaller than the 
width of the ions emissivity curves, or in temperature regions 
where the relative emissivities vary slowly, there is no way 
of constraining the shape of the EMD from measurements 
alone. However, in order to be able to compare different EMD 
solutions, we need to have meaningful confidence intervals. 
Currently, we do not have a theoretical model that describes 
the EMD of a full corona, although ther are theoretical predic- 
tions for the EMD, under various assumptions. We therefore 
make no assumptions as to the shape of the EMD and choose 
to fit a staircase shaped function to allow for local confidence 
intervals estimates. It also makes the equations for the Fe line 
fluxes, which give the best constraints, linear in the fitted pa- 
rameters. 



EMD(T) = {C„ T n <T < T n+ i ; n = 0..JV - 1} 



(5) 



Where N is the number of EM temperature bins, T„ is the lower 
temperature of bin n. C„ > is the parameter to be fitted, which 
is the averaged EMD over the bin. The confidence intervals are 
calculated using the inverse^ 2 distribution, meaning we search 
the parameter space for the x 1 contour that gives a deviation 
from the best-fit that corresponds to the requested confidence 
level. A 1 cr deviation is defined as the contour of y 2 + 1 . The 

^ mm 

selection of the number of bins and their widths is not trivial 
and optimal EMD binning can vary between different spectra. 
The line emissivity curves have considerable widths and some 
extend to temperatures much higher than their peak emissiv- 
ity, resulting in strong negative correlations between the EMD 
in neigbouring bins. Since, as discussed above, we are inter- 
ested in meaningful confidence intervals, we cannot use many 
narrow bins, as this will result in excessive error bars. By mak- 
ing narrow bins they also become closer in temperature, the 
emission contribution of neighbouring bins becomes similar, 
thus increasing the degeneracy. Ultimately, if meaningful con- 
fidence intervals are to be obtained, the number of bins has to 
be kept small, and the tradeoff between temperature resolution 
(number of bins) and constraints (degeneracy of the solution), 
optimized. No smoothing algorithm is applied to the EMD. 

In addition, we fit an EMD parameterized as an exponent 
of a polynomial ensuring that the EMD remains positive: 



EMD(T) = exp 



2 C «Pn(T) 



V n=0 



(6) 



where P„(T) is an n-degree polynomial represented as 
a Chebyshev polynomial for numerical convenience and 
summe d up to N = 7. This is simila r in principle to the method 



used bv lHuenemoerder et alj d2001l) . but with no smoothing. In 
this method, local confidence intervals can not be produced, but 
we use it to verify that our results do not depend on the EMD 
parameterization. 

3.4. Integrated EM 

The physical measurable quantity is the observed line flux, 
which results from the entire plasma (an integration over the 



EMD, eq. Q]). Therefore, the meaningful quantity is the inte- 
gral of the EMD over a range of temperatures. In other words, 
the total EM in that range. Since emissivity curves are smooth, 
uncertainties in the exact way the EM is distributed over a tem- 
perature range, narrower than the emissivity widths, will not 
have a significant effect on the total EM in that range. When 
integrating and taking correlations into account, uncertainties 
caused by the strong correlation between neighboring EMD 
bins that were integrated over, disappear. This results in much 
smaller errors. 

Temperature regions with little emissivity variations lead to 
poor localization of the EMD that in turn create spikes in the 
fitted EMD if no smoothing algorithm is used. Such spikes in 
the solution are likely to appear when too narrow bins are used. 
This reflects the fundamental mathematical instability of the 
solution. The progressively Integrated EM as a function of tem- 
perature (IEM(T)) provides natural smoothing to these spikes 
as the confidence intervals always vary smoothly. When com- 
paring two different EMDs (e.g., flare and quiescence states), 
since we propagate the confidence intervals, the IEM indicates 
with greater certainty where the two solutions differ. It also al- 
lows the comparison of results obtained by different methods 
or by different binning as it clears out correlation related uncer- 
tainties. Over-binning, thus, which can make the EMD mean- 
ingless, will not render the IEM unusable. 

For the staircase EMD, the IEM is tracked at each stage of 
the search for the parameters confidence intervals. This way, 
we get good estimates for the upper and lower limits of the 
IEM curve, inside a^ 2 surface that represents a lcr deviation. 
We also integrate over the EMD from the exponential model to 
verify that the IEM is consistent regardless of the chosen EMD 
parameterization. 

3.5. Abundances 

In order to extract the X/Fe abundance ratios, we calculate the 
non-Fe line fluxes predicted by the fitted, Fe-scaled, EMD. 



pq = A Fe I A z 
fi 4nd 2 \A Fe 



r 

J T 



P ii(T)fq(T)EMD(T)dT 



(7) 



The ratio of the actual measured fluxes to these calculated 
fluxes gives the X/Fe abundance values. Note that for calculat- 
ing abundances, one line from each element is sufficient, while 
for the EMD equations we require two lines of different ion- 
ization degrees. This is why we can estimate abundances for 
elements not used in the EMD fitting. Such is the case with 
Oxygen where we exclude O VII as it forms mostly at tempera- 
tures below our low-T cutt-off, but use O VIII for the abundace. 

Abundance uncertainties are directly related to the specific 
element line flux errors, but they are also indirectly a result of 
EMD uncertainties. We take into account the latter by tracking 
the upper and lower limits on the abundances as we fit the EMD 
and calculate its confidence intervals. This gives the possible 
abundance values within a one cr variation of the EMD. Since 
EMD is determined mostly by Fe and a combination of other 
elements, it is justified to treat the two error contributions as 
uncorrelated. We have also calculated the abundances using the 
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exponential EMD model and verify that the results of the two 
models are in good agreement. 

We do not calculate the absolute abundances, i.e., the abun- 
dances relative to hydrogen. The absolute abundaces are sensi- 
tive to the continuum, which is produced mainly by the less 
co nstrained high-T EMD. This is demonstrated in the work 
of iGu et all J2006) where different instruments observing the 
same source produce different absolute abundances, but similar 
relative abundances. Since we are interested in abundance vari- 
ations during the flare, we use the relative (to Fe) abundances, 
obtained strictly from line fluxes. 



4. Results 

4.1. EMD and integrated EM 

The best-fit EMD for each of the six targets is plotted in fig- 
ure!?] The error bars on the staircase EMD model (eq. (5)) are 
1 cr errors and include correlation uncertainties. The smooth 
dashed line with no confidence intervals is the exponential 
model (eq. (6)). Since the EMD is scaled by the Fe abundance 
and since we do not measure the absolute Fe/H abun dance, we 
assum e in Fig. [4] a solar abundance of 3.24xl(T 5 (Feldman, 
19921) . This has no effect on the following results and is done 
purely for the purpose of setting a reasonable absolute EMD 
scale. 

In all cases, the two models (staircase and exponential) 
agree with each other quite well within the errors over most of 
the temperature range considered. The exponential model tends 
to place a very sharp peak at -8 keV, especially for the flaring 
states. The reason for this is that the constraints on the EMD in 
this region are almost exclusively set by Fe XXIV and Fe XXV. 
In order to set the observed line flux ratio only one temperature 
is needed and the solution tends to a delta function around that 
temperature. When we look at the staircase EMD we can see 
from the large correlation-induced error bars in the last bins, 
that this component can easily be placed in the neighbouring 
bin without changing the resulting spectrum significantly. 

Comparing the flare and quiescence EMD for each target 
we see a general pattern: At temperatures typical of the quies- 
cence state the EMD is similar or increased by some small fac- 
tor and at higher temperatures a new component appears during 
the flare. In CC Eri, the lower- T EMD is increased by a consid- 
erable factor of roughly 3. The HR 1099 flare differs from the 
others by having very little added EMD at temperatures higher 
than those of the quiescence state. 

The integrated EM from zero to kT is plotted in figure[5]for 
all targets with l<x error bars. The smooth dashed line without 
confidence intervals is the integral over the exponential model. 
In the 0.5-2 keV temperature region, where many L-shell emis- 
sivity curves peak, the IEM from both models (staircase and 
exponential) agree very well. At higher temperatures, the lo- 
calization of the EM is not as good and the two models diverge 
slightly, but they always do re-converge at higher temperatures, 
once the integration covers the full range. At 10 keV the total 
EM is remarkably similar for both models and well within the 
error bars. 



In most of the flares (II Peg, OU And, Algol and TZ Crb) 
the IEM of the flare, runs parallel to the IEM of quiescence 
up to ~2 keV in the log-log plot, which indicates that the ther- 
mal structure remains the same or uniformly increased by a 
small factor. In the CC Eri flare, the increase in low-T EM is 
much more significant and reaches a factor of 4.3 at 2 keV, 
where the flare IEM is still increasing while the quiescence 
IEM has nearly reached its maximum value. The HR 1099 flare 
IEM is generally parallel and above the quiescence IEM, except 
around 1.5 keV where they are similar. 

Several theoretical works attempt to predict the EMD of a 
static loop or flaring loops. They commonly predict a power 
law type distribution up to a peak temperature: 



EMD(T) oc T 



(8) 



where a is determined by the details of the cooling processes. 
Note the different definitions of the EMD in various papers. 
For purely radiative cooling, a — —y+ 1 where y is a parameter 
(in the range y « 0+0.5) of the approximated radiative cooling 
function: 



A(T) oc T y 



(9) 



i.e. a = 0.5 - 1.5, while for conduct i ve co oling a « 1.5 and 
for evaporation a ta 0.5 dAntiochosi Il980h . Of special rele- 
vance to this work, where we integrate the spectr a over the en- 
tire flare, is the calculation bv lMewe et al.1 dl997F) . They predict 
for a time-averaged EMD of a quasi-statically cooling flaring 
loop a = 19/8. Since we can resolve the EMD only down to 
0.2 keV, fitting a power law up to the first peak is possible only 
to in those systems that peak is around 2 keV. We also fitted 
the flare in the range of 4-10 keV, where the quiescence emis- 
sion is low. The fits were preformed on the IEM curves as the 
errors there are smaller. The results are shown in table [2] We 
see that a is in the range of 1 .5-2, which is slightly higher than 
expected from radiative and conductive EMD models, signifi- 
cantly higher than expected from evaporation models and sig- 
nificantly lower than expected from the effect of time averaging 
over a single temperature cooling loop. 

4.2. Abundances 

The abundances relative to Fe, as calculated from the stair- 
case model are summarized in Tables lBTHB.6l in the appendix. 
Quoted errors include both flux and EMD uncertainties as ex- 
plained in section 13.51 Abundances measured using the expo- 
nential EMD agree well with these results within a few percent 
of a std-dev. The difference between the EMD representations 
is so small that we only quote the values from the staircase 
model. Since we obtain the element abundance for individual 
ions, elements with more than one ion in the spectrum may 
have several abundance values. The abundance value in the 
table is the statistically weighted average of these individual- 
ion abundances. Note that abundances derived from different 
ions of the same element are expected to be consistent, since 
this was assumed in the EMD reconstruction algorithm (see 
eq. ©). 
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Fig. 4. Emission measure distributions for flare and quiescence states. Solid lines are the binned EMD models with lcr errors. 
Dotted lin es are the expo nent of polynomials EMD models. For scaling purposes only, a solar Fe/H abundance of 3.24xl0~ 5 is 
assumed ( Feldmanl 1992 ). 



5. Discussion 



5.1. Present work 



The EMD plots depict the average variation in the thermal 
structure of the plasma coupled with possible density varia- 



tions. The EMD in each temperature bin represents the aver- 
aged EMD over the bin temperature range. It is apparent that 
the systems that are considered coronally active, namely Algol, 
HR1099 and OU And show a similar EMD shape during qui- 
escence: A gradual increase with kT up to about 2 keV and a 
sharp drop beyond that, becoming negligible at 3-4 keV. This 






kT [keV] kT [keV] 

Fig. 5. The total EM integrated over the EMD from 0.2 keV to kT. Solid line is from the staircase EMD model and includes Icr 
errors. Dotted line is using the exponent of polynomials EMD model. 



is in contrast with the two other systems: TZ CrB and CC Eri, 
which have most of the emission localized around 0.5 keV dur- 
ing quiescence, but still have a small amount of EM up to 2 keV. 
II Peg is a middle case where the EMD peaks around 1 keV and 
a significant portion of the EM is likely to be above 2 keV. 

Looking at figure [5] for most targets, the flare and quies- 
cence IEM curves run nearly parallel (in log scale) to each other 



up to a point where the EM in the flare increases sharply rel- 
ative to quiescence. It is clear that most of the excess EM is 
located at the very high temperatures, typically above 2 keV. 
HR 1099 is the exception to this rule as no significant amount 
of EM is added at high temperatures during the flare. The IEM 
in CC Eri is unique. There is much excess low-T EM (~ factor 
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Table 2. Fitted powerlaws for the EMD. 



System 


quiescence (0.2-2 keV) 
a a red. x 2 


Flare (0.2-2 keV) 
a a red. x 1 


Flare (4-10 keV) 
a cr red. x 2 


II Peg 
OU And 


1.47 0.14 1.33 
1.96 0.24 0.51 


1.48 0.18 0.06 
1.87 0.16 1.4 


0.8" 0.4 0.26 
1.5 0.3 0.2 
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Fig. 6. Flare abundances relative to quiescence abundances as a function of FIR 
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4) in the flare, starting apparently from the lowest temperatures 
observed. 

The similarity in flare and quiescence EMD at temperatures 
below 2 keV, differing only b y a constant factor was also ob- 
served by Nor don et al.l d2006l) in a flare on cr Geminorum. This 
effect was clear from a direct comparison of the flare and qui- 
escence spectra. In that case, the line emission at wavelengths 
above 12 A was increased by a uniform 25%, corresponding 
to a uniform EMD increase by t he same factor under ~2 keV 
temperature. I Audard et al.l (1200 lb also reported a similar effect 
in a flare on HR 1099 detected by XMM-Newton. 

The added EM at low temperatures can be either a real vari- 
ation of the quiescence corona, or an effect caused by the time 
averaging of the flare heating and/or cooling. In other words, 
a narrow temperature component from the flare (moving up or 
down in kT) will mimic a continuous distribution when aver- 
aged over time. Ideally, we would repeat the EMD analysis 
for shorter segments of time to reproduce the time evolution. 
Unfortunately, that also means higher statistical uncertainties 
on the line fluxes, leading to too high uncertainties on the EMD 
to make the comparison feasible. We did extract light curves 
for the brightest lines in the spectrum and looked for a direct 
indication of heating or cooling. No clear conclusion could be 
made. In the case of a narrow temperature component, the cool- 
ing through the low-T range would have to be fast enough in 
order to pass the significant amount of excess EM from high- 
T to below the minimum temperature, without adding much to 
the time averaged EMD. The fact that this added cool EM is 
present even in the II Peg observation, where most of the decay 
phase is not observed, suggests that it is not a time averaging 
effect, but rather a genuine increase in EM over a broad tem- 
perature range in all flares. A note should be added that the 
EMD is scaled by the unknown Fe abundance relative to H. 
Therefore, a possible increase in Fe abundance could explain 
the uniform increase in the cool range EMD. This possibility 
has no effect on, and is consistent with, the FIP related abun- 
dance variations which we discuss below. 

The contrast between the small low-T EM excess and the 
large high-T EM excess could be due to rapid expansion (e.g., 
evaporation, relaxed pinch, loss of magnetic confinement), 
which results in a rapid loss of high-T EM. In this case, the 
flare either reaches low temperatures when its EM is small, or 
the flare originates in a hot environment and therefore cools 
down only to the (high) ambient temperature. 

Table [3] summarizes the values of the total integrated EM 
(last data point in figure[5]). We see that in most cases the flares 
cause an increase in EM by a factor of 2-4. In HR 1099, inter- 
estingly, it increased by only 25%, while in CC Eri it increased 
by a factor of -25, although the total EM of CC Eri is still 
very small relative to the other systems. This means that the 
quiescent corona of CC Eri is very small and the emission in 
the flare could be attributed entirely to the flare heated plasma. 
In the other systems we observe a mix of flare plasma with a 
background of the strong quiescent corona. 

Overall, there are no dramatic variations in abundances, but 
several elements do show a tendency of decreased abundances 
relative to Fe during the flares, mainly O and Ne . A similar ef- 
fect has also been observed in a flare on cr Gem dNordon et al 



2006). It should be noted that in those flares where the low-T 
EM (kT < 2 keV) during the flare is not significantly larger than 
in quiescence, the measured flare abundances are practically a 
weighted average of the abundances in the flaring plasma and 
the background quiescence plasma. 

Figure [6] shows the flare to quiescence abundance ratios 
plotted as a function of FIR Since the abundances are plotted 
relative to Fe it means that if the Fe/H abundance has changed 
in the flare, all the ratios in the figure will be multiplied by a 
uniform factor of Fe qu i es /Fefi are . This may re-scale the plot, 
but will not change the pattern or our conclusions. We see that 
in all flares, with the exception of TZ CrB, high FIP elements 
show reduced abundances. For the low FIP elements the pic- 
ture is not as clear as the elements with FIP lower than Fe: Al 
and Ca, have large error bars. Still, elements with FIP of 8 eV 
and below seem to vary together with Fe (flare/quies is ~1). 
This behaviour is similar to the solar FIP effect. TZ CrB is the 
exception, with no apparent abundance effects as all abundance 
ratios are consistent with unity. 

This is somewhat surprising as other works ha ve shown that 



increa sed coronal activity leads to the IFIP effect. lAudard et al 



(2003) have examined abundances on RS CVn systems and 
found that high FIP elements were increasingly over-abundant 
as the typ ical temperature (use d as an indication of activity) 
increased. iTelleschi et al.l (120051) have found a similar effect in 
solar like stars where the FIP effect switched to IFIP as activity 
increased. For stars with an IFIP corona, chromospheric com- 
position would appear to be FIP in comparison. The change 
in abundances during the flare indicates that the excess plasma 
is not heated coronal plasma, but more likely heated chromo- 
spheric plasma. This would mean that, if flares are at all respon- 
sible for the IFIP effect, the selective element transport would 
have to operate during the cooling or post-flare stage since we 
detect a low-FIP enriched composition during the flare itself. 

From looking at the light curves in figure Q] we clearly see 
two types of flare behaviours: Symmetric flares with a sharp 
peak and a rapid decay (CC Eri & Tz CrB), and asymmetric 
flares that rise fast, but decay slowly, which creates a broad 
peak (II Peg, HR1099, OU And). Algol, at first glance, appears 
to belong to the first kind, but the flare occurre d just as the 
sy stem was coming out of eclipse. As noted by IChung et el 



sy 
(2i 



2004), it is likely that if most of the emission originates from 
near chromospheric level, the rise phase in the X-ray light curve 
is governed by the exposure of the flaring region behind the 
eclipsing companion and not by the rise of the flare. Both of 
the rapid decay flares occurred on the cooler dwarf stars, while 
the long duration flares are all on the giant or sub-giant stars. 
However, CC Eri flare shows the clearest FIP effect in the flare 
relative to quiescence, similar to the long duration flares, while 
TZ CrB shows no abundance variations. The sample is too 
small to draw any substantiated conclusions, but it seems that 
there is no clear correlation between the type of flare and abun- 
dance effects. 
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Table 3. Total integrated EM" from 0.2 to 10 keV and averaged 
luminosities* in the 1.8-20 A range. 





Quiescence 




Flare 


EM Ratio 


Flare Int. 


Total Net Flare 




EM 


u c 


EM 


U c 


Flare/Quies. 


Time 


Energy 


Traget 


[10 53 cm" 3 ] 


[10 30 erg s~'] 


[10 53 cm" 


3 ] [10 30 erg s" 1 ] 




[ks] 


[10 34 erg] 


II Peg rf 


1 47 +0 - 32 


9.8+0.2 


■x 7+0.9 
-0.35 


22.5+0.2 


2.5 


22.7 


28.7±0.6 


OU And 


15 25 +4J 


40.6+0.8 


63.2?™ 


95.5+1.3 


4.1 


39.5 


216±6 


Algol 


3 14+ 037 

->- 1 ^-0.43 


7.69±0.06 


8-3!!;? 


12.8+0.2 


2.6 


13.7 


7.0+0.2 


HR 1099 


7 4+0.5 
^ -0.3 


9.61±0.08 


•'•<•::: 


10.7±0.1 


1.25 


15.9 


1.85+0.2 


TZCrB 


1 1Q+0.07 
1 • 1 7 -0.03 


2.92±0.02 


5 r - 7 

J ' 1 -0.6 


9.5+0.2 


4.3 


6.2 


4.1+0.1 


CCEri 


"•""-0.007 


0.223±0.003 


2 +0 - 2 

^•"-0.2 


2.42±0.04 


25.3 


8.0 


1.76±0.03 


" Using the staircase EMD 



Luminosities are averaged over the time intervals as indicated on the 
light curves 

c Luminosity of the averaged spectra during quies./flare states 
d Observation ended while excess flux was still high. 



5.2. Comparison With Previous Works 

Several of the observations used in this work (table [TJ were 
published before by other researchers, although for the most 
part the focus was not on the flare effects. Also, the methods 
were different than ours and the presentation of results makes 
comparisons difficult. We bring here a short summary of the 
results of those individual works as well as references to other 
relevant high-resolution X-ray observatio ns of these targets. 
Obs. 1451 of II Peg was analyzed by IHuenemoerder et al 



(2001) who found that during the flare, at temperatures of 
log T — 7.3 to log T = 8.0 K, a large EM component was added 
while the cooler emission was hardly changed. We find a sim- 
ilar pattern, but the cooler emission tends to be slightly higher 
during the flare, which could be due to Fe abundance variation 
as previously discussed. The abundance ratio of Ne/Fe relative 
to solar drops from 22+6 preflare to 17+5 (a factor of 0.77) 
during the flare, which is only slightly more than the present 
ratio of 0.86+0.05. 

Obs. 1892 of OU And has not been fully analyzed to date. 
The system has also been ob served by XMM -Newton during a 
quiescent state as reported bv lGondoin (|2003j). 



Obs. 604 of Algol was analyzed by IChung et elj ([2004). 
They investigated slight line shifts and showed that the X-ray 
emission is dominated by the secondary star in the system. The 
corona is likely to be asymmetric and located closer toward the 
center of mass of the system. There is no special treatment of 
the flare. Algol is a well studie d target in quiesce nce by the 
high resolution instruments. lSchmitt & Nessl d2004l) used a line 
fitting method on a different Chandra LETG observation, to 
get a smooth EMD solution for Algol in its quiescence state. 
Depending on the details of the fitting, they get various EMD 
solutions with peaks at 1-2 keV, consistent with ours. A flare 
on Algol has also been observed during an eclipse by XMM- 
New ton and allow e d fo r spatial information to be extracted, 
see: ISchmitt et alj (12003b . Only limited spectral information 



was used. 

Obs. 62538 of HR 1099 was performed as part of a multi- 
wavelength campaign involving Chandra, EUVE, HST and 
the VLA. Results from this collaboration were published by 
Avres et al.1 (1200 lb . however no EMD or abundances were pub- 



lished. HR 1099 was also ob served during a larg e flare with 
XMM-Newton as reported by lAudard et al. ( 200 lb who found 
an abundance enhancement of the low-FIP elements Fe and Si 
during the flare, while high-FIP Ne remained constant. (The 
quiescent abundances were reported to have an overall IFIP 
trend.) This FIP bias during the flare is similar to what we find 
for HR 1099 and in genera l for five of the six flares we have an- 
alyzed. |Osten_etalJ ([2004J) used EUVE observations and found 



an HR 1099 quiescence EMD in agreement with the one pre- 
sented here, showing a broadly distributed EM in the range of 
0.4 keV to 2.5 keV where their plot is cut-off. 

Obs. 15 of TZ CrB which was part of a campaign that 
included Chandra, EUVE, and the VLA, was analyzed by 



Osten et al.l Q2003). They report an overall increase in abso- 
lute abundances and no FIP related pattern during either qui- 
escence or flares. In our work absolute abundances (relative 
to H) were not calculated, but we also find no FIP pattern in 
the flare relati ve to quiescence. W e also detect a Ni XIX line at 
12.43 A, while lOsten et alJd2003h claim the absence of Ni XIX 
lines. Though this Ni line is somewhat blended with Fe lines, 
the combined emissivity of Fe contributes less than 50% of 
the total flux in the feature according t o our flux measure ment 
method. An EMD was constructed by lOsten et al.l ([2003), but 
no error bars are given, making a comparison difficult. Their 
EMD, similar to ours, features two major components, one at 
0.5 keV which exists also in quiescence and the other at about 
3 keV which emerges with the fl are. This t a rget h as also been 
observed with XMM-Newton by ISuh et al.l d2005l) . They con- 
firm the absence of a clear FIP bias. 

Obs. 6132 and 4513 of CC Eri were not published yet. CC 
Eri was previously observed in low resol ution by AS CA, si- 
multaneously with EUVE and the VLA dOsten et al. 
No flares were observed. 



6. Conclusions 

We have analyzed six large X-ray flares observed with Chandra 
on six different systems. Emission measure distribution and in- 
tegrated EM were calculated, including well-localized confi- 
dence intervals allowing for an unambiguous comparison be- 



2002). 
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tween flare and quiescence states. Relative abundances were 
measured from the EMD and compared between states. We 
verified that the derived abundances do not depend on the rep- 
resentation of the EMD. Integration over both staircase and 
exponential EMDs produce the same total EM to high accu- 
racy, showing that the inevitable degeneracy of the EMD so- 
lution has little effect on the total EM. Thus, we conclude that 
in contrast with the local EMD(T), the EM integrated from a 
low temperature to a (continuously varying) high temperature 
is a convenient quantity useful for comparison between differ- 
ent emission states and different targets. 

During the six flares analyzed in this work, the EMD at 
temperatures below ~2keV appears to be similar to that dur- 
ing quiescence with a small, roughly uniform enhancement ob- 
served during the flare. In five of the six flares, the added EM 
is predominantly at temperatures of kT > 2 keV. The total EM 
is increased by a factor of 2-4 for most flares, but by a factor of 
-25 for the CC Eri flare and only by 25% for HR 1099. Five out 
of the six targets show a statistically significant flare FIP bias in 
which the high-FIP to low-FIP abundance ratios decrease dur- 
ing the flares. The exception is the TZ CrB flare that showed no 
statistically significant abundance variations. We conclude that 
in our sample, flaring activity tends to evaporate plasma with 
abundances biased toward a (solar) FIP effect and not an IFIP 
effect. Note that this conclusion does not require the knowledge 
of and therefore is independent of photospheric or quiescent 
coronal abundances. 

The different element composition observed during the 
flare and the significant increase in total EM indicate that the 
flaring plasma is likely to be heated chromospheric plasma, 
rather than locally heated coronal plasma. This might support 
the chromospheric evaporation scenario. On the other hand, our 
results are inconsistent with previous suggestions that flaring 
activity and chromospheric evaporation produce an IFIP effect. 
If flares are indeed responsible for the transition from FIP to 
IFIP coronal composition in very active coronae, and in order 
to be consistent with our observations, the fractionation mech- 
anism must operate during the cooling stage or after the flare. 
If what we observe in flares is indeed photospheric composi- 
tion plasma heated to coronal temperatures, flares may provide 
a mean to measure photospheric abundances in active stars. On 
the other hand, it could be the flare itself that introduces the FIP 
bias by selectively injecting low FIP material into the corona. 
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Table A.l. Measured line fluxes in 10 4 photons s _1 crrT 2 of II Peg and OU And 



Wave 



Quiescence 
Flux Err 



II Peg 

Flare 
Flux Err 



Flare/Quies. 
Ratio Err 



Quiescence 
Flux Err 



OUAnd 

Flare 
Flux Err 



Flare/Quies. 
Ratio Err 



15.01 
14.21 
13.51 
12.84 
12.28 
11.77 
11.00 
10.64 
1.85 
1.78 
21.60 
18.97 
13.45 
12.13 
9.17 
8.42 
7.76 
7.17 
6.65 
6.18 
5.04 
4.73 
3.95 
3.73 
3.18 
3.02 
12.43 



1.18 

0.64 

0.54 

0.47 

0.69 

0.50 

0.42 

0.50 

0.04 

0.00 

3.26 

16.75 

3.40 

10.73 

0.42 

0.76 

0.04 

0.09 

0.53 

0.57 

0.31 

0.45 

0.140 

0.061 

0.000 

0.000 

0.031 



0.17 

0.15 

0.14 

0.13 

0.11 

0.09 

0.07 

0.07 

0.14 

0.27 

2.20 

1.00 

0.23 

0.32 

0.05 

0.07 

0.04 

0.04 

0.06 

0.07 

0.10 

0.12 

0.073 

0.076 

0.060 

0.079 

0.073 



1.45 

0.83 

0.72 

0.68 

0.80 

0.66 

0.59 

1.06 

0.38 

0.15 

2.36 

20.32 

3.96 

13.83 

0.48 

1.23 

0.00 

0.21 

0.60 

0.99 

0.35 

0.63 

0.130 

0.132 

0.131 

0.073 

0.091 



0.17 
0.16 
0.15 
0.15 
0.12 
0.10 
0.07 
0.09 
0.17 
0.28 
2.03 
1.02 
0.23 
0.33 
0.06 
0.08 
0.04 
0.06 
0.06 
0.08 
0.12 
0.14 
0.087 
0.098 
0.082 
0.110 
0.086 



1.23 
1.30 
1.35 
1.43 
1.16 
1.31 
1.40 
2.12 
10.13 

0.72 
1.21 
1.17 
1.29 
1.16 
1.62 
0.00 
2.33 
1.13 
1.73 
1.14 
1.40 
0.93 
2.17 



0.22 
0.39 
0.45 
0.49 
0.26 
0.29 
0.28 
0.35 
37.70 

0.79 
0.09 
0.10 
0.05 
0.21 
0.18 
1.13 
1.25 
0.17 
0.25 
0.54 
0.47 
0.79 
3.16 



2.98 7.66 



0.48 
0.27 
0.26 
0.27 
0.36 
0.23 
0.32 
0.50 
0.16 
0.03 
0.05 
1.57 
0.17 
0.96 
0.13 
0.40 
0.007 
0.029 
0.23 
0.35 
0.068 
0.13 
0.046 
0.011 
0.017 
0.019 
0.000 



0.08 
0.06 
0.06 
0.06 
0.05 
0.03 
0.03 
0.03 
0.07 
0.13 
0.77 
0.25 
0.05 
0.07 
0.02 
0.03 
0.014 
0.016 
0.02 
0.03 
0.036 
0.04 
0.028 
0.031 
0.027 
0.031 
0.028 



1.00 

0.43 

0.40 

0.40 

0.61 

0.43 

0.55 

1.18 

0.83 

0.30 

0.00 

2.31 

0.16 

1.16 

0.22 

0.80 

0.01 

0.10 

0.32 

0.69 

0.19 

0.25 

0.053 

0.081 

0.061 

0.117 

0.045 



0.13 
0.10 
0.09 
0.08 
0.08 
0.05 
0.05 
0.06 
0.13 
0.21 
1.11 
0.40 
0.07 
0.10 
0.03 
0.05 
0.02 
0.03 
0.03 
0.05 
0.06 
0.07 
0.047 
0.055 
0.046 
0.056 
0.046 



2.09 
1.63 
1.55 
1.47 
1.70 
1.88 
1.71 
2.34 
5.09 
9.58 
0.01 
1.47 
0.94 
1.21 
1.73 
2.02 
1.45 
3.38 
1.42 
1.97 
2.85 
1.95 
1.15 
7.10 
3.67 
6.20 



0.43 
0.52 
0.50 
0.43 
0.31 
0.36 
0.21 
0.19 
2.21 
40.17 
20.65 
0.35 
0.50 
0.14 
0.37 
0.22 
4.27 
2.13 
0.20 
0.21 
1.77 
0.86 
1.25 
19.59 
6.53 
10.63 



NOTE: Fluxes of H-like ion lines include both transitions of the unresolved Ly-a doublet. 
Fluxes of He-like ion lines include only the resonant transition. 

Fluxes of L-shell Fe ion lines include all transitions within ±0.03 A of the specified wavelength. 
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Table A.2. Measured line fluxes in 10" 4 photons s" 1 crrT 2 of Algol and HR 1099 



Algol 







Quiescence 


Flare 


Flare/Quies. 


Quiescence 


Flare 


Flare/Quies. 


Ion 


Wave 


Flux 


Err 


Flux 


Err 


Ratio 


Err 


Flux 


Err 


Flux 


Err 


Ratio 


Err 


Fe XVII 


15.01 


4.77 


0.18 


5.55 


0.35 


1.17 


0.09 


4.12 


0.22 


5.73 


0.31 


1.39 


0.11 


Fe XVIII 


14.21 


2.67 


0.15 


2.73 


0.28 


1.02 


0.12 


2.09 


0.19 


2.49 


0.29 


1.19 


0.17 


FeXIX 


13.51 


2.33 


0.15 


2.50 


0.28 


1.07 


0.14 


1.70 


0.16 


1.93 


0.22 


1.14 


0.17 


FeXX 


12.84 


2.12 


0.14 


2.39 


0.27 


1.13 


0.15 


1.78 


0.22 


1.71 


0.21 


0.96 


0.17 


FeXXI 


12.28 


2.25 


0.13 


2.73 


0.24 


1.21 


0.13 


1.94 


0.14 


1.95 


0.19 


1.00 


0.12 


Fe XXII 


11.77 


1.82 


0.09 


2.04 


0.17 


1.12 


0.11 


1.18 


0.10 


1.35 


0.14 


1.14 


0.15 


Fe XXIII 


11.00 


1.79 


0.08 


2.67 


0.15 


1.50 


0.11 


1.23 


0.08 


1.64 


0.11 


1.33 


0.12 


Fe XXIV 


10.64 


2.52 


0.08 


4.76 


0.17 


1.89 


0.09 


1.35 


0.08 


2.00 


0.12 


1.49 


0.13 


FeXXV 


1.85 


0.43 


0.12 


2.02 


0.35 


4.69 


1.58 


0.22 


0.12 


0.27 


0.20 


1.22 


1.12 


Fe XXVI 


1.78 


0.05 


0.17 


0.67 


0.46 


12.94 


44.59 


0.00 


0.19 


0.06 


0.32 


- 


- 


O VII 


21.60 


4.25 


1.68 


0.02 


2.88 


0.00 


0.68 


1.48 


1.79 


0.82 


2.60 


0.55 


1.88 


O VIII 


18.97 


11.73 


0.57 


11.88 


1.15 


1.01 


0.11 


26.27 


0.97 


24.48 


1.22 


0.93 


0.06 


NelX 


13.45 


2.12 


0.15 


1.85 


0.26 


0.87 


0.14 


4.39 


0.21 


4.54 


0.29 


1.04 


0.08 


NeX 


12.13 


7.70 


0.20 


8.56 


0.39 


1.11 


0.06 


16.04 


0.33 


15.86 


0.43 


0.99 


0.03 


Mg XI 


9.17 


0.79 


0.06 


0.92 


0.10 


1.16 


0.15 


0.74 


0.06 


0.86 


0.09 


1.18 


0.15 


Mg XII 


8.42 


1.82 


0.07 


2.88 


0.14 


1.58 


0.10 


1.46 


0.08 


1.64 


0.12 


1.12 


0.10 


Al XII 


7.76 


0.07 


0.03 


0.15 


0.07 


2.10 


1.42 


0.23 


0.04 


0.16 


0.06 


0.69 


0.29 


Al XIII 


7.17 


0.25 


0.04 


0.35 


0.08 


1.45 


0.42 


0.19 


0.05 


0.25 


0.07 


1.31 


0.47 


Si XIII 


6.65 


0.99 


0.06 


1.42 


0.11 


1.43 


0.14 


0.83 


0.06 


0.89 


0.08 


1.07 


0.13 


Si XIV 


6.18 


1.51 


0.07 


2.65 


0.15 


1.75 


0.13 


0.95 


0.07 


1.32 


0.11 


1.38 


0.15 


SXV 


5.04 


0.60 


0.09 


0.79 


0.19 


1.32 


0.37 


0.39 


0.10 


0.75 


0.16 


1.90 


0.63 


SXVI 


4.73 


0.52 


0.09 


0.87 


0.20 


1.66 


0.48 


0.53 


0.11 


0.66 


0.17 


1.25 


0.41 


Ar XVII 


3.95 


0.21 


0.06 


0.23 


0.13 


1.11 


0.72 


0.31 


0.08 


0.33 


0.11 


1.06 


0.44 


Ar XVIII 


3.73 


0.10 


0.06 


0.31 


0.15 


3.13 


2.55 


0.20 


0.08 


0.18 


0.12 


0.87 


0.67 


CaXIX 


3.18 


0.21 


0.06 


0.38 


0.12 


1.84 


0.77 


0.19 


0.07 


0.24 


0.10 


1.21 


0.68 


CaXX 


3.02 


0.04 


0.06 


0.28 


0.16 


6.83 


11.42 


0.20 


0.09 


0.00 


0.11 


0.00 


0.53 


NiXIX 


12.43 


0.37 


0.073 


0.52 


0.15 


1.41 


0.49 


0.32 


0.08 


0.23 


0.12 


0.74 


0.42 



HR 1099 



NOTE: Fluxes of H-like ion lines include both transitions of the unresolved Ly-a doublet. 
Fluxes of He-like ion lines include only the resonant transition. 

Fluxes of L-shell Fe ion lines include all transitions within ±0.03 A of the specified wavelength. 
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Table A.3. Measured line fluxes in 1(T 4 photons s _1 cirT 2 of TZ CrB and CC Eri 



Wave 



Quiescence 
Flux Err 



TZ CrB 

Flare 
Flux Err 



Flare/Quies. 
Ratio Err 



CC Eri 

Quiescence Flare Flare/Quies. 

Flux Err Flux Err Ratio Err 



15.01 
14.21 
13.51 
12.84 
12.28 
11.77 
11.00 
10.64 
1.85 
1.78 
21.60 
18.97 
13.45 
12.13 
9.17 
8.42 
7.76 
7.17 
6.65 
6.18 
5.04 
4.73 
3.95 
3.73 
3.18 
3.02 
12.43 



12.76 

5.04 

3.65 

2.78 

2.47 

1.43 

0.97 

0.87 

0.06 

0.01 

2.48 

19.92 

4.25 

9.28 

1.54 

1.67 

0.14 

0.10 

1.16 

0.78 

0.38 

0.14 

0.11 

0.02 

0.07 

0.01 

0.65 



0.18 
0.12 
0.10 
0.09 
0.08 
0.05 
0.04 
0.04 
0.04 
0.07 
0.35 
0.57 
0.12 
0.15 
0.04 
0.04 
0.02 
0.02 
0.03 
0.03 
0.04 
0.04 
0.02 
0.02 
0.02 
0.02 
0.04 



14.20 
6.91 
4.29 
4.12 
3.49 
1.97 
2.13 
4.93 
1.68 
0.17 
0.71 
32.78 
5.38 
15.54 
1.70 
4.07 
0.17 
0.44 
1.84 
2.54 
1.28 
1.07 
0.27 
0.36 
0.40 
0.35 
1.10 



0.83 
0.69 
0.54 
0.51 
0.44 
0.29 
0.24 
0.30 
0.57 
0.84 
6.56 
2.95 
0.61 
0.78 
0.22 
0.25 
0.12 
0.15 
0.19 
0.24 
0.35 
0.37 
0.24 
0.26 
0.22 
0.27 
0.30 



1.11 

1.37 
1.18 
1.48 
1.41 
1.37 
2.19 
5.68 

27.42 

12.13 
0.29 
1.65 
1.26 
1.67 
1.10 
2.43 
1.17 
4.52 
1.58 
3.28 
3.38 
7.79 
2.41 

18.00 
5.83 

25.25 
1.69 



0.07 
0.14 
0.15 
0.19 
0.18 
0.21 
0.26 
0.43 
20.26 
81.06 
2.65 
0.16 
0.15 
0.09 
0.15 
0.16 
0.88 
1.74 
0.17 
0.35 
0.99 
3.39 
2.21 
24.46 
3.59 
42.87 
0.47 



2.49 
0.86 
0.57 
0.51 
0.41 
0.26 
0.23 
0.22 
0.03 
0.01 
3.66 
14.75 
2.96 
4.96 
0.27 
0.28 
0.02 
0.02 
0.37 
0.28 
0.18 
0.10 
0.04 
0.02 
0.00 
0.01 
0.16 



0.09 
0.06 
0.05 
0.04 
0.03 
0.02 
0.02 
0.02 
0.03 
0.05 
0.41 
0.51 
0.08 
0.10 
0.01 
0.02 
0.01 
0.01 
0.02 
0.02 
0.02 
0.02 
0.01 
0.01 
0.01 
0.01 
0.02 



7.72 
3.43 
2.79 
2.21 
2.03 
2.53 
3.44 
6.98 
3.15 
0.67 
7.28 
40.15 
5.10 
16.81 
1.01 
3.37 
0.07 
0.30 
1.93 
4.83 
2.05 
2.53 
0.70 
0.40 
0.33 
0.22 
0.60 



0.75 
0.64 
0.48 
0.42 
0.40 
0.28 
0.24 
0.31 
0.59 
0.77 
7.85 
3.62 
0.55 
0.74 
0.15 
0.25 
0.10 
0.12 
0.16 
0.26 
0.33 
0.37 
0.22 
0.23 
0.20 
0.22 
0.25 



3.09 
3.97 
4.91 
4.35 
4.95 
9.66 
15.16 
32.05 



1.99 
2.72 

I. 72 
3.39 
3.71 

12.12 
3.66 

13.52 
5.18 

17.55 

II. 31 
24.91 
18.08 
20.63 

20.97 
3.87 



0.32 
0.80 
0.93 
0.90 
1.06 
1.35 
1.59 
3.14 



2.16 
0.26 
0.19 
0.16 
0.58 
1.21 
5.25 
7.37 
0.49 
1.41 
2.34 
6.53 
8.32 
18.37 

33.24 
1.73 



NOTE: Fluxes of H-like ion lines include both transitions of the unresolved Ly-a doublet. 
Fluxes of He-like ion lines include only the resonant transition. 

Fluxes of L-shell Fe ion lines include all transitions within ±0.03 A of the specified wavelength. 
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Appendix B: Measured Abundances 



Table B.l. II Peg abundances relative to Fe during flare and 
quiescence. 





Quies. 


Flare 


Ratio 




El. 


X/Fe 


Error 


X/Fe 


Error 


Flare/Quies. 


Error 


O 


133 


10 


97.2 


9.6 


0.73 


0.09 


Ne 


51.56 


2.5 


44.5 


1.6 


0.864 


0.053 


Mg 


2.92 


0.26 


2.75 


0.19 


0.94 


0.10 


Al 


0.30 


0.13 


0.332 


0.01 


1.10 


0.57 


Si 


2.39 


0.22 


2.015 


0.15 


0.84 


0.10 


S 


1.60 


0.36 


1.11 


0.21 


0.70 


0.21 


Ar 


0.57 


0.28 


0.31 


0.16 


0.54 


0.39 


Ca 


N/A 


N/A 


0.28 


0.17 


N/A 


N/A 


Ni 


0.02 


0.05 


0.048 


0.045 


2.4 


6 



Table B.2. OU And abundances relative to Fe during flare and 
quiescence. 





Quies. 


Flare 


Ratio 




El. 


X/Fe 


Error 


X/Fe 


Error 


Flare/Quies. 


Error 


O 


22.4 


4.3 


16.1 


3.0 


0.72 


0.19 


Ne 


7.47 


0.62 


4.00 


0.36 


0.536 


0.065 


Mg 


1.98 


0.16 


1.80 


0.13 


0.906 


0.097 


Al 


0.118 


0.064 


0.169 


0.051 


1.44 


0.89 


Si 


1.59 


0.12 


1.244 


0.086 


0.781 


0.079 


S 


0.46 


0.14 


0.41 


0.090 


0.88 


0.33 


Ar 


0.16 


0.11 


0.118 


0.064 


0.74 


0.65 


Ca 


0.10 


0.12 


0.149 


0.062 


1.4 


1.8 


Ni 


N/A 


N/A 


0.038 


0.040 


N/A 


N/A 



Table B.3. Algol abundances relative to Fe during flare and 
quiescence. 



Quies. Flare Ratio 

El. X/Fe Error X/Fe Error Flare/Quies. Error 



o 


22.0 


1.8 


16.62 


2.2 


0.75 


0.12 


Ne 


8.90 


0.25 


6.77 


0.32 


0.760 


0.042 


Mg 


1.507 


0.055 


1.513 


0.073 


1.004 


0.061 


Al 


0.168 


0.027 


0.167 


0.036 


0.99 


0.27 


Si 


1.171 


0.046 


1.198 


0.062 


1.023 


0.066 


S 


0.500 


0.06 


0.400 


0.068 


0.80 


0.16 


Ar 


0.172 


0.048 


0.133 


0.050 


0.77 


0.36 


Ca 


0.215 


0.06 


0.184 


0.054 


0.85 


0.34 


Ni 


0.063 


0.013 


0.080 


0.023 


1.27 


0.45 



Table B.4. HR 1099 abundances relative to Fe during flare and 
quiescence. 





Quies. 


Flare 


Ratio 




El. 


X/Fe 


Error 


X/Fe 


Error 


Flare/Quies. 


Error 


O 


74.1 


3.4 


51.0 


3.1 


0.69 


0.05 


Ne 


25.23 


0.67 


19.0 


0.6 


0.75 


0.03 


Mg 


1.77 


0.09 


1.58 


0.10 


0.89 


0.07 


Al 


0.29 


0.05 


0.24 


0.05 


0.82 


0.22 


Si 


1.26 


0.07 


1.20 


0.08 


0.95 


0.08 


S 


0.65 


0.11 


0.75 


0.12 


1.15 


0.27 


Ar 


0.48 


0.10 


0.35 


0.11 


0.72 


0.27 


Ca 


0.44 


0.14 


0.31 


0.14 


0.68 


0.38 


Ni 


0.065 


0.017 


0.04 


0.02 


0.61 


0.35 



Table B.5. TZ CrB abundances relative to Fe during flare and 
quiescence. 





Quies. 


Flare 


Ratio 




El. 


X/Fe 


Error 


X/Fe 


Error 


Flare/Quies. 


Error 


O 


25.1 


1.7 


26.0 


3.6 


1.04 


0.16 


Ne 


7.95 


0.38 


8.15 


0.44 


1.025 


0.073 


Mg 


1.66 


0.07 


1.71 


0.10 


1.026 


0.076 


Al 


0.138 


0.015 


0.172 


0.054 


1.25 


0.41 


Si 


1.25 


0.06 


1.176 


0.090 


0.943 


0.085 


S 


0.492 


0.055 


0.53 


0.12 


1.08 


0.28 


Ar 


0.240 


0.055 


0.142 


0.080 


0.59 


0.36 


Ca 


0.29 


0.09 


0.177 


0.080 


0.60 


0.33 


Ni 


0.057 


0.005 


0.076 


0.021 


1.33 


0.38 



Table B.6. CC Eri abundances relative to Fe during flare and 
quiescence. 





Quies. 


Flare 


Ratio 




El. 


X/Fe 


Error 


X/Fe 


Error 


Flare/Quies. 


Error 


O 


72.3 


3.3 


34.6 


4.1 


0.48 


0.06 


Ne 


21.2 


0.5 


9.8 


0.5 


0.465 


0.025 


Mg 


1.49 


0.07 


1.34 


0.10 


0.90 


0.08 


Al 


0.13 


0.04 


0.10 


0.04 


0.78 


0.37 


Si 


2.17 


0.09 


1.49 


0.08 


0.689 


0.046 


S 


1.26 


0.15 


0.76 


0.09 


0.60 


0.10 


Ar 


0.40 


0.12 


0.17 


0.05 


0.44 


0.19 


Ca 


0.08 


0.14 


0.10 


0.05 


1.2 


2.1 


Ni 


0.074 


0.012 


0.077 


0.033 


1.05 


0.48 



